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ABSTRACT 

Near-infrared spectroscopy using the new Palomar Integral Field Spectrograph in- 
dicates the presence of an AGN (active galactic nucleus) in the ultraluminous infrared 
galaxy IRAS 08311—2459. The high-velocity wings of the Paschen-a hydrogen recom- 
bination line are seen to be spatially unresolved, and with no positional offset between 
red and blue high velocity emission. The [Si VI] coronal line, with a 167 eV excita- 
tion potential, is also spatially unresolved with a velocity width comparable to that of 
the broad component of the Paa emission. The low velocity component of the Paa 
emission is seen to be rotating, and is extended over ~ 2 kpc. Molecular hydrogen is 
also extended, and elongated along the maximum velocity gradient, sharing the same 
rotation profile as the narrow Paa emission. The simple picture in agreement with the 
observations is that of an AGN surrounded by a rotating disk of star formation. These 
observations lend strength to the usage of [Si VI] as a diagnostic of AGN activity in ul- 
traluminous infrared galaxies, and also highlight the utility of integral field spectroscopy 
in elucidating subtle morphological differences in line emitting regions. 

Subject headings: galaxies: active — galaxies: individual (IRAS 08311—2459) — galaxies: 
infrared — galaxies: starburst 

1. Introduction 

Ultraluminous infrared galaxies (ULIRGs) are among the most luminous sources in the uni- 
verse, with infrared luminosities of L{ r > 1O 12 L0. The bulk of the energy emitted from these 
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sources emerges in the far-infrared, suggesting that the light we see is thermally processed by warm 
dust in the galaxies. As such, the optically thick dust prohibits direct viewing of the mechanism 
responsible for the extreme power generation. Quasars (or more generally QSOs — quasi stellar 
objects) and massive starbursts are perhaps the only phenomena capable of producing the high 
luminosities observed in these systems. Recent spectroscopic programs at infrared wavelengths 
have shown that the great majority of ULIRGs with Lj r < 2 x 1O 12 L are powered predominantly 
by starbursts rather than by AGN (Genzel et al. 1998; Rigopoulou et al. 1999; Murphy et al. 1999; 
Veilleux, Sanders, & Kim 1997, 1999). At luminosities higher than this, the AGN fraction rises to 
between one-third and one-half. 

IRAS 08311-2459 is a ULIRG at cz = 30150 km s -1 (1 arcsec = 1.7 kpc) with L ir = 2.5 x 
1O 12 L , and is classified as a "warm" ULIRG by the definition of Sanders et al. (1988). This galaxy 
is selected from the 2Jy ULIRG sample of Strauss et al. (1990, 1992), and is a single-nucleus galaxy 
surrounded by tidal debris (Murphy et al. 1996). The near-infrared spectrum of IRAS 08311—2459 
(Murphy et al. 1999) shows emission line features suggestive of AGN activity, though a scenario 
involving strong shocks and high velocity outflows could also explain these spectral features. This 
galaxy was selected for the present study because it is the only galaxy in a volume-limited sample 
of 33 ULIRGs (Murphy et al. 2000) showing strong [Si VI] emission, and one of two displaying 
any obvious high-velocity Paa emission. Calculated quantities involving luminosities and physical 
scales assume a cosmology with Hq = 75 km s _1 Mpc -1 and go = throughout this paper. 

Near-infrared integral field spectroscopy, effectively providing simultaneous velocity-resolved 
imaging capabilities of multiple emission lines, allows subtle morphological differences between 
various line emitting regions to be elucidated. In this paper we present integral field observations 
of the Paa, H 2 1-0 S(3), and 167eV excitation [Si VI] emission lines in IRAS 08311-2459, which 
together form a consistent picture of a galaxy containing a central AGN surrounded by a starburst 
disk. 



2. Observations and Data Reduction 

Observations were made using the Palomar Integral Field Spectrograph (PIFS) situated at 
the //70 focus of the 200-inch Hale Telescope. This instrument produces simultaneous spectra 
for eight slits in a contiguous two-dimensional field of view. A description the PIFS instrument, 
along with general observing and data reduction procedures can be found in Murphy, Matthews, 
& Soifer (1999). IRAS 08311-2459 was observed on the night of 25 March 1999. The R « 1300 
resolution mode (Av ~ 225 km s -1 ) was used to obtain two spectra centered on the redshifted 
Paa, and H2 1-0 S(3)+[Si VI] lines, respectively. Separate sky exposures were alternated with the 
on-source integrations, with integration times of 300 s. A positional dither pattern was employed 
for the sequence of integrations enabling recovery of seeing-limited spatial resolution in the cross- 
slit direction. The pixel scale for these observations is 0.167 arcsec pixel -1 . An offset field star 
was used for auto guiding, with image motion compensated by driving the active secondary mirror. 
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Guided observations of a nearby star, with 10 s exposures, accompanied the spectral observations 
for the purpose of evaluating the point spread function (PSF) . Wavelength calibration is provided 
through a combination of OH airglow lines (Oliva & Origlia 1992) and arc lamp spectra taken at 
the time of observation. All wavelengths are referred in air units. Atmospheric opacity and spectral 
flat-fielding are compensated simultaneously using the light from HR 3862, a 4.9 mag GO V star, 
spread uniformly across the field of view. At the time of observation, the G star was at an airmass 
of ~ 1.95, well matched to the airmass of the science observations. 

For these observations, the 5 / /4x9 / /6 field of view was oriented with the long axis at a position 
angle of 90°. The Paa line was observed with 1800 s of on-source integration time at an average 
airmass of 2.00. The H2+[Si VI] lines were observed for 2100 s of on-source time at an average 
airmass of 1.96. The seeing for the Paa dataset, as determined via imaging of the PSF star, was 
measured at l'/OO midway through the observation, and 1'/15 at the end. The telescope was then 
focused, resulting in a PSF full-width at half-maximum (FWHM) of l'.'05. Observations of the 
H2+[Si VI] line immediately followed, with the PSF measuring 0?85 midway through, and / /95 at 
the end of the observation. The measures of seeing should not be taken entirely at face value, as 
variability over the time scales of the spectral observations limits the accuracy of the very short 
PSF exposures at estimating the average seeing conditons. Furthermore, the sliced PSF images are 
adequate for assessing the approximate size of the seeing disk, but of insufficient quality to permit 
characterization of subtle distortions such as may be caused by slight telescopic astigmatism. All 
observations presented here were made in clear photometric conditions. Photometric calibration 
was performed via K s imaging of the galaxy and of the faint standard star 9143 from Persson et 
al. (1998), using a flip-in mirror in front of the grating, thereby establishing the continuum flux 
density at 2.155 fim. The continuum slopes observed in the Paa and H2+[Si VI] spectra were used 
to estimate the continuum level for the Paa spectrum. 

Data reduction consists of first performing the sky subtraction, interpolating static bad pix- 
els and cosmic ray artifacts, division by the G star spectrum, and multiplication by a blackbody 
spectrum matched to the G star's temperature of 5930 K. The H2+[Si VI] spectrum was multiplied 
by a template G3 V spectrum from Kleinmann Sz Hall (1986), smoothed and resampled to the 
PIFS resolution, in order to remove the Br7 stellar absorption line introduced in the division by 
the HR 3862 spectrum. No calibrator absorption lines of any significance effect the Paa spectrum. 
Spatial and spectral distortions are corrected using previously generated distortion maps appropri- 
ate for the particular grating setting. Co-registration of the eight slits in the spatial dimension is 
based upon observation of the G star with its light extended perpendicular to the slit pattern by 
chopping the telescope secondary mirror in a triangle-wave pattern. The two-dimensional spectra 
from the eight slits are placed into the three-dimensional datacube according to the positional dither 
pattern, with a common wavelength axis established by the calibration lines. Residual OH airglow 
lines are removed by subtracting a scaled version of the raw sky spectrum, with typical scalings 
of ~2% in absolute value, sometimes as large as 5%. Photometric variability among individual 
integrations is compensated by small scaling adjustments such that the object flux is consistent 



-4 - 



from one integration to the next. A more complete description of the data reduction procedures 
may be found in Murphy, Matthews, & Soifer (1999). 

In the data presented here, the pure continuum images are constructed directly from the 
integral field datacube, and constructed from spectral regions free of line emission as well as OH 
airglow emission. For each spatial pixel, a linear fit is made to the line- free portion of the continuum 
spectrum. The line images are formed by subtracting this continuum fit from the datacube, then 
summing over a small range in the spectral dimension. In this way, the line images simulate 
continuum-free narrow band images at arbitrary wavelengths — or velocities — within the datacube. 
The co-registration of the continuum and line images is implicit, as the images are derived from 
the same dataset. 

3. Results 

3.1. Paa Emission Line 

Figure 1 shows the appearance of the Paa spectrum (A res t = 1.87510/xm), taken from a 
/ /83xO / . / 83 square aperture centered on the continuum peak of the galaxy. The resolved veloc- 
ity profile of the line shows broad wings around a moderately narrow peak. The wings appear 
asymmetric, with more emission blue-ward of the line center. Part of this asymmetry likely comes 
from the He I line at a rest wavelength of 1.8689 / um, with an estimated contribution of ~ 5% 
of the Paa flux, as is characteristic in ULIRGs (Murphy et al. 1999). Subtracting the expected 
He I contribution with the same velocity profile as Paa almost completely removes the blue excess, 
resulting in a Paa profile consistent with being symmetric. A multiple Gaussian line fit, whose com- 
ponents are displayed in Figure 1, fits a well-centered 800 km s _1 broad component plus a narrow, 
160 km s _1 wide component to the Paa line, with comparable total fluxes. See Appendix A and 
Table 1 for details on the line fits. A small amount of excess emission appears in the fit residuals 
of Figure 1 on the blue side of the Paa line. This emission only amounts to ~4% of the flux of the 
broad Gaussian component, which was not substantial enough to encourage a broad component fit 
that had any asymmetric quality (i.e., blue-shifted line center). Therefore, while this emission is 
real, it is of little importance in comparison to the bulk of the broad emission. Throughout this 
paper we refer to this 800 km s -1 emission feature as the "broad" component, though it does not 
represent a canonical broad line region. We simply use this term to indicate velocities larger than 
would be expected in ordinary galaxy kinematics, such as rotation. 

The continuum morphology is presented in Figure 2, along with images of three different spec- 
tral ranges within the Paa line profile. Referring to the line peak at 2.06345/im as the zero velocity 
reference, the three continuum-subtracted narrow band images correspond to central velocities of 
—425, 0, and +425 km s -1 , each 250 km s _1 wide. The narrow line image is more compact than the 
continuum, with characteristic FWHM values of 1'.'55 and l'/75, respectively. The line emission is 
also slightly elongated compared to the round continuum, with an ellipticity (major-to-minor axis 
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Fig. 1. — Nuclear spectrum of the Paa line extracted from a 07 83x0? 83 aperture in the datacube. 
Positions of the He I and H2 1-0 S(4) lines are marked at the systemic velocity determined from 
Paa. In addition to the spectral data, the continuum fit and Gaussian fit components are shown, as 
is the residual to the functional fits. The broad component of the Paa emission is clearly seen, and 
is rather symmetric after accounting for the He I line. The absolute flux density scale corresponds 
to light within the aperture only. Only 18% of the total continuum flux from a 5" x 5" box centered 
on the galaxy falls within this aperture. Aperture corrections for the line emission can be found in 
Table 1 of the Appendix. 
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Fig. 2. — Morphologies of the continuum and of the Paa emission in three velocity bins, centered at 
v = and v = ±425 km s" 1 , each 250 km s" 1 wide. All show rather round, symmetric morphologies, 
though the broad emission is significantly more compact than either the narrow Paa or continuum 
emission. The FWHM of the broad line emission measures compared to 1//55 and l'/75 for 

the narrow Paa and continuum emission, respectively. North is up, and east is to the left in all 
images. Contours are placed at y/2 multiplicative intervals. 

ratio) of ~ 1.1. The two images of the line wings are significantly more compact than either the 
continuum or narrow line Paa images, and are in fact consistent with the estimated l'.'l seeing limit 
at the time of observation. No significant difference in shape, size, or center can be distinguished 
between the red and blue line wing images, though they each exhibit an ellipticity of ~ 1.12. The 
observed ellipticity could be the result of a slightly defocused, astigmatic telescope. Such a distor- 
tion would not noticeably distort the much larger continuum shape, and likewise would not fully 
account for the similar elongation of the narrow Paa image. 

With confidence in the generic morphological character of the line emission, we may turn to a 
more detailed representation of the line extents by plotting the line emission FWHM as a function 
of wavelength. Figure 3 shows the measured spatial FWHM in both the x and y dimensions of 
the array, corresponding to the east-west and north-south directions. Each data point represents 
a measurement in the spatial plane at the location of a single spectral pixel. Because a spectral 
resolution element spans roughly four pixels, adjacent points in Figure 3 are not independent. The 
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Fig. 3. — Spatial line extent as a function of wavelength for the Paa line, showing the spatially 
resolved nature of the narrow line emission in contrast to the spatially unresolved high velocity 
regions. Two plots are shown, representing measures in the orthogonal array directions. The size 
of the continuum source is shown at the top of the plot. The estimated range in seeing is indicated, 
as measured from PSF images intermixed with the spectral data. The profile at bottom shows the 
continuum subtracted spectrum of the Paa line from Figure 1 for immediate comparison to the 
spatial extents. The increase in spatial width at the left edge of the plot is attributable to the 
He I line. Line centers are marked for the systemic velocity, and a velocity scale is provided for 
convenient reference. 

horizontal dotted lines delimit the range of atmospheric seeing measured in conjunction with the 
spectral observations, with a dashed line indicating the mean of these two values. The extended 
nature of the narrow line emission stands out in clear contrast to the seeing-limited extent of the 
line wings. The increase in width at the left edge of the plot coincides with the red side of the He I 
emission line. Due to the low signal-to-noise ratio in the region of the He I line, it was not possible 
to extend the plot data to shorter wavelengths. The continuum FWHM, measured beyond the line 
wings, is represented by the dashed line at top, and the continuum-subtracted line profile is plotted 
for convenient reference at bottom. 



3.2. H 2 10 S(3) and [Si VI] Emission Lines 



The integral field spectrum of IRAS 08311-2459 at 2.15/im contains the Br<5, H 2 1-0 S(3), 
and [Si VI] emission lines at rest wavelengths of 1.94456/xm, 1.95702/xm, and 1.96287/xm, respec- 
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Fig. 4. — Nuclear spectrum of the H 2 1-0 S(3) + [Si VI] line region, extracted from a / . / 67x0 / /67 
aperture in the datacube. The Bt5 line is also seen. In addition to the spectral data, the continuum 
fit and Gaussian fit components are shown, as is the residual to the functional fit. The [Si VI] line 
is clearly broader than the other lines in the spectrum, and comparable in strength to the H2 line. 
Line centers are marked at the systemic velocity. The absolute flux density scale corresponds to 
light within the aperture only. Only 13.5% of the total continuum flux from a 5" x 5" box centered 
on the galaxy falls within this aperture. Aperture corrections for the various emission lines can be 
found in Table 1 of the Appendix. 

tively. Figure 4 shows the one-dimensional spectrum, extracted from a / . / 67x0 / . / 67 square aperture 
centered on the continuum peak. The Br<5 line is largely ignored here, as the Paa line is a far more 
effective probe of hydrogen recombination, and the combination of Br<5's weakness plus its spectral 
proximity to Paa limits its usefulness as a measure of extinction. The 1'.'35 spatial extent of the 
Br 5 line is perfectly consistent with the Paa extent, accounting for the seeing difference between 
the observations. The [Fe II] line at 1.9670 ^m, seen in some ULIRGs (Murphy et al. 1999), does 
not appear in the spectrum, though the Br7 feature from the atmospheric calibrator coincides with 
this same wavelength, making the presence of [Fe II] uncertain. 

All lines in the H2+[Si VI] spectrum are spectrally resolved, with the H2 line resembling the 
atomic hydrogen lines in terms of observed line width. Though the H2 line is blended with the 
[Si VI] line, it is known to be symmetric by comparing with the H2 1-0 S(l) line profile (Murphy 
et al. 1999). The [Si VI] line is significantly broader than its neighboring lines, with a FWHM of 
~ 1000 km s -1 . See Appendix A and Table 1 for more information on the functional fits to these 
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lines. 

Figure 5 displays the continuum and narrow band line images of IRAS 08311—2459 in H2 1-0 
S(3) and [Si VI] emission. The continuum morphology appears different in Figures 2 and 5, mainly 
due to 20% better seeing and telescope focus in the latter (see Section 2) . In Figure 5 the continuum 
appears to have a weak extension to the SW. The line images are clearly more compact than the 
continuum, and the near perfect symmetry of the [Si VI] emission suggests that the structure seen 
in the continuum is real, and simply not resolved in Figure 2. Indeed the central contours of the 
continuum image in Figure 5 appear much more symmetric than the outer contours. The line 
images are created in the same manner as those in Figure 2, with the spectral range encompassing 
points within ±200 km s~ x of the line center. The [Si VI] emission line is slightly more compact 
than the neighboring H2 line. An elliptical Gaussian fit to each of these line images yields FWHM 
sizes (average of major and minor axes) of l'.'l for H2 and 0f!9 for [Si VI], with ellipticities (major- 
to-minor axis ratios) of 1.30 and 1.03, respectively. The H2 elongation appears to be real, with a 
position angle on the sky of ~ 30° — similar to the disk orientation implied by the rotation axis, 
discussed in Section 3.3. 

A more detailed look at the difference in spatial distribution between H2 and [Si VI] can be 
seen in Figure 6, which, like Figure 3, plots the spatial FWHM of the line distribution in the east- 
west and north-south directions as a function of wavelength. Referring to the line profile at the 
bottom of Figure 6, it is clear that throughout the wavelength range corresponding to the broad 
[Si VI] emission, both the x and y widths are consistent with the seeing limit, indicating that the 
entire [Si VI] emission line is spatially unresolved. The similarity in widths in the two orthogonal 
directions further demonstrates the symmetry of the PSF during the observation. On the short 
wavelength side of the [Si VI] emission, the spatial profile broadens, marking the location of the 
H2 line. Here, the north-south extent is larger than the east-west extent, owing to the elongated 
nature of the emission at a 30° position angle. 

3.3. Rotation of H2 & Narrow Paa 

Taking advantage of the coexistence of two-dimensional spatial plus spectral information, we 
can construct velocity fields of the line emitting gas. At each spatial pixel, the wavelength of the 
peak line emission is computed, and converted to a velocity. Figure 7 displays the velocity fields 
obtained in this manner for IRAS 08311—2459 in both the Paa and H 2 1-0 S(3) lines. Both lines 
exhibit a velocity structure indicative of pure rotation. The color scale is the same for both maps, 
through which it is seen that the rotation is not only aligned along the same axis for the two lines, 
but with roughly the same amplitude. The rotation axis is at a position angle of ~ 110°, implying 
a projected disk major axis of ~ 20°, consistent with the observed elongation observed for H2. 
Figure 7 also shows the rotation curves obtained for both lines, in a / .'67 wide aperture through 
the galaxy center at a 20° position angle. Though the molecular hydrogen data has much more 
scatter than does the Paa curve, it is clear that the two constituents — the narrow Paa and the H2 
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Fig. 5. — Continuum and line morphologies for the H2 and [Si VI] lines. The line images are each 
400 km s _1 wide, centered on the line peaks. The [Si VI] emission is completely unresolved, as well 
as symmetric. The H2 line is slightly resolved, and extended along a 30° position angle, roughly 
consistent with the disk direction inferred from the rotation. The continuum shape differs slightly 
from that seen in Figure 2, in that a small extension appears to the SW. The seeing conditions were 
better for this image, and the compact symmetry of the [Si VI] line suggests that this extension is 
real. Orientation and contours are as described for Figure 2. 




Fig. 6. — Spatial line extent as a function of wavelength for the H2 1-0 S(3) and [Si VI] lines, 
showing the spatially unresolved nature of the [Si VI] emission in contrast to the spatially resolved 
H2 emission. See the caption for Figure 3 for a description of plot features. The spectrum at 
bottom is a continuum subtracted version of that seen in Figure 4. 



-12- 



200 

100 





>, 



o -100 

73 

> 



-200 



_>> 

'u 
_o 

CD 
> 



u 
cu 
w 
u 







■ 


^^^^ 



200 






2 


100 






xsec 


1 




100 
200 


1 




-1 
-2 



■2-10 1 2 
arcsec 



V. 



hi 



-2-10 1 2 
arcsec 



E 



> 



00 
50 




- + Paa 
A H 



8 -50 



-1 00 
-1 50 



-NNE+ + a 



A 



A 



A 
A 



A 





A 



A,+ 



A 



SSW : 



-2-10 1 2 
displacement (") 



Fig. 7. — Velocity fields of both the Paa and H2 lines. To the left are color images representing 
the velocity of the peak emission, with red indicating redshift, etc. A clear rotational signature 
is seen in both. To the right is a position-velocity plot for both species along a C'67 "slit" at a 
position angle of 20°, corresponding to the angle of maximum velocity gradient. The Paa points 
(plus symbols) trace out a very clean rotation pattern, and the noisier H2 points (triangles) follow 
the same profile. 



emission — share the same physical rotation. 



4. Discussion 



Combining visible light spectroscopic surveys of ULIRGs (e.g., Kim, Veilleux, & Sanders 1998; 
Veilleux, Kim, &; Sanders 1999) with infrared observations (Veilleux, Sanders, & Kim 1997, 1999), 
it is found that 20-25% of ULIRGs can be characterized as containing AGN. The visible studies 
generally utilize the [O IIl]/H/3 line ratio as the primary diagnostic in differentiating Seyfert (AGN) 
galaxies from starburst-dominated galaxies, while infrared identification of AGN are generally based 
on the presence of broad atomic recombination lines, and occasionally of [Si VI]. ULIRGs are known 
to sometimes harbor galactic-scale superwinds (Heckman, Armus, & Miley 1990), which could 
potentially contribute to the observed broad line wings and to [Si VI] emission through powerful 
shock excitation. The present integral field data for IRAS 08311—2459 argue for a simple picture 
of a compact AGN surrounded by a more diffuse starburst — at least some of which is occurring in 
a rotating disk around the nucleus. While massive outflows and expanding shells are present in 
some ULIRGs, they may not be ubiquitous. The PIFS data on IRAS 08311—2459 highlights the 
important role that integral field spectroscopy can play in untangling the dynamics and excitation 
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of the circumnuclear gas in ULIRGs. 

Without considering the contribution of the He I line to the Paa line profile, one might conclude 
that the apparently blue-asymmetric profile of Paa indicates the presence of outflow phenomena, 
as is often the case for asymmetric line profiles (e.g., Heckman, Armus, & Miley 1990; Wang et al. 
1997). Knowledge of the He I line on the blue side of Paa, with expectations on typical ULIRG line 
strengths (Murphy et al. 1999, 2000) allow us to understand the Paa profile as the combination of 
a symmetric Paa line mixed with a similarly shaped He I line. The Gaussian line fits give the He I 
line 4% of the total Paa flux — consistent with expectations from the ULIRG population (Murphy 
et al.). Perhaps the strongest evidence that the He I line is contributing to the total line profile 
comes from the observation of increased spatial extent coincident with the red side of the He I line, 
shown in Figure 3. 

While we can not rule out the presence of a very compact nuclear outflow, the symmetry of the 
Paa line profile resulting from proper accounting of the He I flux, together with the fact that both 
sides of the profile's broad wings are spatially unresolved (< 500 pc) in the full two-dimensional 
sense, argue that the broad emission does not arise from an outflow, or wind, phenomenon. Ad- 
ditionally, the red and blue wing components show no relative displacement greater than the Cf.'05 
level, or ~ 100 pc, further supporting this conclusion. If a wind is to account for the high-velocity 
Paa emission, the lack of significant extinction to the red emission (Ay < 0.5 mag based on line 
symmetry) must be explained, as must the fact that the physical scale is smaller than the size of 
the narrow line region, which is spatially resolved. We believe the broad emission arises from the 
near vicinity of an obscured AGN, and that this component is separate from the rotating, narrow 
emission discussed at the end of this section. 

The presence of [Si VI] in the spectrum of IRAS 08311—2459 is a clear indicator of high 
energy processes at work in the galaxy. [Si VI] is not commonly seen in ULIRGs, even those 
exhibiting broad Paa emission (Murphy et al. 1999). Veilleux, Sanders, & Kim (1997, 1999) see 
evidence for [Si VI] in a few ULIRGs, though blending with the comparably strong H2 1-0 S(3) 
line prohibits accurate estimates of flux contributions. The [Si VI] emission in IRAS 08311—2459 
measures 6 x 10 Lq in luminosity, uncorrected for extinction. This is among the most luminous 
[Si VI] emission line yet observed, although relative to the bolometric luminosity only slightly more 
luminous than that found in nearby Seyfert galaxies (Marconi et al. 1996). The [Si Vl]-to-bolometric 
luminosity ratio for quasars is the same as that seen in IRAS 08311—2459, based on the template 
QSO spectrum in Murphy et al. 

Main sequence stars are incapable of producing significant flux at 167 eV — the requirement 
for creating Si 5+ . Only ionizing radiation from an AGN or fast shocks can account for the [Si VI] 
emission (Contini & Viegas 1992). It is difficult to differentiate between these excitation processes 
observationally, though a diagnostic utilizing the [Fe VII] A0.6087 / um line was explored by Marconi 
et al. (1996), in which all instances of [Si VI] detections in Seyfert galaxies are attributed to 
photoionization. 
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The integral field data may contribute a clue to the nature of the [Si VI] excitation. The 
molecular hydrogen emission — often associated with collisional excitation — has a significantly dif- 
ferent morphological character than the [Si VI] emission. Namely, the H2 line originates from a 
region roughly one kpc larger than that of [Si VI]. The underlying morphology of violent shocks in 
major mergers is not clearly understood. However, if the [Si VI] line was shock excited — perhaps 
tracing out a physical boundary along a shock front — then any molecular hydrogen excited by this 
same shock phenomenon would arise from the near vicinity of the shock front (Draine & McKee 
1993), thereby sharing similar spatial structure with the [Si VI] emission. Naturally, the very strong 
shocks capable of producing Si +5 would probably dissociate H2 molecules, but the precursor and 
post-shock mechanisms may be capable of producing significant H2 emission. In order for both the 
[Si VI] and H2 emission to be shock-excited requires two separate shock components with vastly 
different velocities and spatial scales. While we can in no way rule out this possibility with the 
current data, a simpler explanation has the [Si VI] emission originating from photoionization by a 
central AGN, with the H2 emission excited by other means. 

Further differentiating the H2 distribution from that of [Si VI] is the fact that the H2 emission 
is extended along the direction of maximum velocity gradient — the inferred disk axis. A simple 
quadrature decomposition of the H2 l'.'Ox 1'.'3 FWHM ellipse with a Of! 9 seeing disk yields an intrinsic 
aspect ratio of greater than 2:1, with a long axis extent of about 1.6 kpc. The rotation curve derived 
from the H2 emission is in complete agreement with that seen for Paa (Figure 7) , suggesting that the 
narrow Paa and H2 emission arise from the same rotating disk. Similar observations of elongated 
molecular gas concentrations coincident with disk orientation and rotation are also seen in CO data 
of luminous infrared galaxies (Bryant & Scoville 1999). 

A picture emerges wherein both the H I and H2 narrow line emission arises from a rotating, 
star forming disk. Alternatives for production of the narrow line emission include X-ray photoion- 
ization of the disk from the central source, or photoionization of the gas from "precursor" radiation 
emanating from a powerful shock front. Both of these alternatives suffer from the same difficulty of 
requiring the ionizing radiation to penetrate almost 1 kpc through the disk in order to produce the 
large scale emission observed. If indeed star formation is occurring on a scale of ~ 1.5 kpc — as the 
H2 emission indicates — then IRAS 08311—2459 boasts star formation on a significantly larger scale 
than the 100-300 pc typically seen among ULIRGs in mid-infrared studies (Soifer et al. 2000). 

If the narrow, rotating component of the Paa line is produced by young stars in a disk, the 
implied star formation rate (SFR) is ~ 75 Mq yr _1 , assuming an intrinsic Ha/Paa line flux ratio of 
8.6 (Osterbrock 1989, Case B with n e = 10 cm -3 and T = 10000 K) and a conversion between Ha 
luminosity and SFR following Kennicutt (1983). This estimate only includes the flux in the narrow 
Gaussian component to the Paa line, as portrayed in Figure 1. Of course, this conversion takes 
into account neither the extinction at Paa, nor the contribution of the central AGN to the narrow 
Paa line flux — both of which may be significant in this system, yet would affect our SFR estimate 
in opposite senses. Unpublished visible light spectra of IRAS 08311—2459 indicate an extinction to 
the disk of Ay k, 1.1 mag, comparing the Paa and H/3 emission, leading to 0.16 mag of extinction 
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at Paa. The AGN contribution may exceed this value, such that the above SFR estimate likely 
represents an upper limit. 

If the entire far-infrared luminosity in IRAS 08311—2459 were attributable to star formation, 
then one may compute the expected total star formation rate by two methods. The first, from 
Scoville & Young (1983), puts a lower bound on the bolometric luminosity by assuming that the 
luminosity is dominated by stars converting protons to He via the CNO cycle, and that the rate 
of mass reduction is the rate at which stars are forming. This approach yields an expected SFR of 
~ 200 M Q yr _1 . The other method, from Hunter et al. (1986), integrates the total stellar luminosity 
output to produce an estimated SFR of 650 M yr _1 , assuming all of the stellar luminosity is 
reprocessed into infrared emission {(3 = 1 in their model). Therefore, it is likely that either there 
is a large amount of star formation that remains obscured at 2 /xm, or that the hidden AGN is 
a substantial contributor to the infrared luminosity in this system. Given the presence of the 
unresolved, relatively strong [Si VI] emission we favor the latter interpretation, though with an 
estimated SFR of ~75 M yr -1 , star formation could clearly contribute a significant portion of the 
total luminosity. 

5. Conclusions 

Integral field spectroscopy of the suspected AGN-powered ULIRG, IRAS 08311—2459, provides 
important constraints on the nature of the line emitting regions in this galaxy. The evidence 
presented can be reasonably explained by the presence of a bona-fide quasar embedded within a 
disk of star formation in the post-merger system. Both components contribute significantly to the 
total energy production in the galaxy. The salient points derived from these data are: 

1. The low velocity emission of Paa, along with H2, has been identified with a rotating disk 
of material with a diameter of ~ 1.5 kpc. It is likely that widespread star formation is 
responsible for the observed narrow line emission from this region. 

2. Moderately broad emission is seen in both Paa and [Si VI], with a projected FWHM« 1000 
km s -1 . The broad emission is spatially unresolved (< 500 pc) in both species. Both lines 
are spectrally symmetric, and there is no observed spatial offset between the red and blue 
wings of the broad Paa emission. Outflow phenomena are therefore unlikely to account for 
the high velocity emission, with an AGN being the preferred source of the [Si VI] and broad 
Paa emission. 
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A. Functional Line Fits 

Functional line fits have been applied to the extracted nuclear spectra of IRAS 08311—2459, 
and displayed in Figures 1 and 4. In each case, the continuum is fit with a linear function, and 
all lines except for Paa are fit with a single Gaussian profile. Attempting to fit Paa with a 
single Gaussian resulted in an obviously inadequate fit, but the resulting parameters are useful in 
fitting weaker atomic lines. The signal-to-noise ratio in the weaker lines does not justify the use 
of multiple component fits, yet their underlying profiles are expected to closely resemble Paa. The 
single-fit Paa has a central velocity of 30147 km s _1 and a FWHM of 415 km s _1 ; accounting for 
an instrumental resolution of 230 km s _1 , this translates to an intrinsic line width of 345 km s _1 . 
These numbers become the constraints used, where necessary, for the low ionization atomic lines. 

The Paa spectrum contains three observable emission features, fit with four Gaussian profiles — 
two for Paa itself. It was found to be necessary to fix the He I line center and width according to 
the single-Gaussian Paa values as discussed above, leaving only the flux as a free parameter. The 
two components in the Paa line fit were unconstrained, and tended to repeatably settle to the same 
parameters despite factors of two modifications to the initial estimated values. The H2 1-0 S(4) 
line was constrained to share the 380 km s _1 line width found for the H2 1-0 S(3) line, as discussed 
below. The plot of fit residuals at the bottom of Figure 1 shows that the two-Gaussian fit for Paa 
does a passable job, but excess emission on the blue side well above the noise level demonstrates 
the approximate nature of the fitting business. For this reason, the separate components in the 
Paa fit are not to be taken as individual, distinct physical components, but rather as a tool for 
gaining a qualitative feel for the nature and symmetry of the broad emission. Table 1 lists the line 
properties and their fit parameters. 

The H2+[Si VI] spectrum contains three lines, each fit with a single Gaussian profile. The fit 
was performed placing no constraints on the three parameters of each line fit. Figure 4 shows the 
result of the combined fit, and Table 1 gives the properties of the individual line fits. 
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Table 1. Fit Parameters and Line Fluxes 



Line 


cz - 30150 
(km s- 1 ) 


FWHM a 
(km s- 1 ) 


Equiv. Width 
(nm) 


Aperture b 
Correction 


Total Flux c 
(xl0~ 18 W m- 2 ) 


He I 


d 


339 d 


0.78 


0.21° 


3.6 


Paa narrow 


-9 


156 


9.3 


0.19 f 


45.9 


broad 


16 


779 


12.0 


0.30 f 


34.1 


H 2 1-0 S(4) 


7 


376 d 


0.48 


0.26 c 


1.7 


BrS 


-14 


297 


0.97 


0.15 f 


3.5 


H 2 1-0 S(3) 


28 


384 


2.6 


0.21 f 


7.4 


[Si VI] 


-66 


1066 


5.6 


0.30 f 


10.8 



a Deconvolved via simple quadrature by the instrument resolution of 230 km s 1 in the Paa 
spectrum, and 220 km s _1 in the H2 + [Si VI] spectrum. 

b Ratio of flux within extraction aperture (i.e., spectra presented in Figures 1 and 4) to total 
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spatial characteristics (as determined from Figures 3 and 6), this ratio is measured directly from 
the datacube. 
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